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Abstract. Most of the solar physicists use very expensive software for data 
reduction and visualization. We present hereafter a reliable freeware solution 
based on the Python language. This is made possible by the association of 
the latter with a small set of additional libraries developed in the scientific 
community. It provides then a very powerful and economical alternative to 
other interactive data languages. Although it can also be used for any kind of 
post-processing of data, we demonstrate the capabities of such a set of freeware 
tools using THeMIS observations of the second solar spectrum. 



1. Introduction 



The Python language ( van Rossum |[l990l ) is increasingly popular in the scientific 



and, more specifically, in the astronomical community. For instance, it has 
been already adopted by very large communities such as the Space Telescope 
Science Institute (STScI) or by the ALMA project. It is also very strong among 
the various tools used and developed in the frame of the Virtual Observatory 
endeavour (see |http : //www . ivoa . net/| for instance). 

In combination with a small number of specific libraries, we explain briefly 
here how it can provide a very powerful interactive data language which can be 
run on a number of operating system platforms, such as Linux, Mac OS X or 
Windows. The set of Python libraries we shall describe herafter does not present 
an interest limited to data reduction since it can also be perfectly used for any 
kind of post-processing like, for instance, the one of output from numerical 
simulations. 

For the spectropolarimetric data reduction we performed here, we used a 
combination of tools from the libraries PyFITS, numarray (now NumPY/SciPy) 
and matplotlib. 

Using data collected by G. Molodij and F. P aletou in May 2000 at the 
THeMIS solar telescope (iPaletou fc Molodii \\2m\ ). we demonstrate the capa- 



bilities of such tools by extracting the so-called 2nd spectrum of the Sri spec- 
tral line at 460.7 n m close to the solar limb, as an illustrative example (e.g., 
Stenfio et al.~lll997l ^. 



By reprocessing such data with our Python-based tools, we cou ld easily 
reproduce the results published elsewhere ( Trujillo Bueno et al. |[200ll ). 
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Figure 1. a: — A image of the fractional linear polarization Q(A) //(A) in the 
spectral domain around 460.7 nm. 



2. Necessary resources 

In order to deal with FITS format files, the STScI developed - and maintains 
- the PyFIT^ library. It is very easy to install and to use since it is very well 
documented. It allows for both the reading of such files, entirely or by slices, 
and for the generation of new FITS files; working with headers is also very easy. 
However its use implies that another library such as numarray or NumPy (see 
below) for handling multi-dimensional arrays is already available. 

Indeed, for our purpose vector calculations with multi-dimensional arrays 
have been performed with the numarray library. With the latter comes along 
also a number of high-level numerical tools allowing for linear algebra, statisti- 
cal analysis, fast fourier transforms, convolutions or interpolations for instance. 
However, since numarray will not be supported anymore after 2007, we wish 
to warn the reader to use inst ead, from now on, the NumPj|3 package for such 
scientific calculations (see also Oliphanti 12006 !). An on-line cookbook and very 



useful documentation can be found at jhttp : / /www . sc ipy . org/, 

Finally, for graphical output and figures saving, we used the matplotlit0 
library for 2D plots. The default GUI is very convenient, allowing for a posteriori 
interactive work on the image, such as area selection and zooming. Unlike other 
software, it is also very easy with matplotlib to export images into the most 
useful formats. The q uality of the outp ut is perfectly suitable for publication 
(see e.g., the figures in lLeger et al.l [20071 ). 



^http: //www, stsci . edu/resources/sof tware.heLrdweire/pyf its] 

^http : / / nvunpy . scipy . org/ 1 

• ^http : //matplotlib . sourcef orge . net/] 
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Figure 2. Mean linear polarization profile in the 460.7 nm spectral domain. 
The profile results from the averaging along the slit of the map displayed in 
Fig. 1. The depolarization of a nearby Fei line at pixel « 33 can also be 
noticed. 

The freeware set we adopted was installed and used by us without any 
difficulties together with several Linux distributions as well as Mac OS X. 



3. Capabilities and results 

The raw data consisted in a time-sequence of 200 frames taken with the slit par- 
allel to the solar limb at fi ~ 0.1 while modulating in time through a sequence of 
4 independent polarization states in order to perform full-Stokes measurements. 

Demodulation involving inverse or pseudo-inverse matrix calculation was 
very easily coded with the linear algebra package of NumPy. Our older routines 
(written in IDL) for aligning spectral lines were also very quickly re-coded us- 
ing convolution with a finite width kernel for bisector search, and cubic spline 
or Fourier interpolations for shifting the profiles line by line. The adequate 
functions are in the f ftpack and signal packages of SciPy. 

In Fig. 1, we plot the (x — A) image corresponding to the fractional linear 
polarization Q(X) /I(X) obtained by recomb ining the two polarized beams of the 
MTR@THeMIS flPaletou Molodii l[200ll ') taken at 460.7 nm. The strong peak 
in the Sri spectral line corresponds to ~ 1.5%. Fig. 2 makes this reading easier 
since it is the mean profile obtained by averaging the polarized signal over the 
rows of the previous image. 
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It is in agreement with previous values obtained at very high polarimet- 
ric sensitivity with the Z IMPOL I polar imeter attached to the NSO/Kitt Peak 
MacMath-Pierce facihty (jStenflo et ah IIT997i ). 

As with other data language, these packages permit to build one's own 
collection of specific functions and, therefore makes it possible for any user to 
constitute its own library. The use of these public and private resources can be 
made both from scripts and/or interactively, using a command line. 

One of the advantages of using Python is flexibility. Several other high-level 
plotting libraries exist such as DISLIJN0, or VTK0 for 3D graphics, for instance. 
And should one be unsatisfied of some of the functions of matplotlib, it would 
be very easy to use instead a number of functions from such other libraries. 



4. Conclusions 

The Python language and the numerous scientific and graphic libraries which are 
being developed for, already provide very valuable and powerful tools for data 
analysis in astrophysics. It is well supported by increasingly larger communities 
so, shifting to such freeware tools appear to us as both reasonable and economical 
(not to say legal too...) options yet. 

Our experience with the reduction of THeMIS spectropolarimetric data 
since 1999 made rather fast, on the timescale of a few weeks only, the conversion 
of our former software written in IDL to the above-mentioned Python-based 
resources. We could reprocess without any difficulties whatsoever old data which 
lead to published results, and we now reduce our new data with those numerical 
tools. And finally, we definitely adopt a proselytizing attitude in favour of them. 

Acknowledgments. Rafael Manso is highly acknowledged for his constant 
interest during this quest for freeware solutions. 
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